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The inclusive and exclusive B decays into the rf meson plus others are investigated 
in a model based on the QCD anomaly. The invariant mass distribution is discussed. 
The constraint of the effective coupling is obtained from the data of the exclusive 
decay modes. The branching ratio of J/ip — > rfr}^ is predicted to be as large as 
5.4 x 10 -5 , which can be tested in the forthcoming CLEO-c experiments. 

1. Introduction. 

A few years ago, unexpected large branching ratios of B decaying into final states with an rf 
meson such as B — > rj'X s and B — > rfK were observed by the CLEO collaboration [1, 2] and 
recently confirmed by BaBar[3] and Belle [4]. This stimulated many theoretical activities in 
understanding the special role of the rf meson in B decays. As the contribution of traditional 
four quark operators from the effective Hamiltonian in the standard model (SM) is far below 
the data[5, 6], various exotic mechanisms were introduced such as a large coupling between 
the gluon and rf through the QCD anomaly[7-ll], intrinsic cc content inside 77' [12, 13] 
and positive interference between several contributions in the SM[14, 15] et.al. The large 
contribution arising from new physics was also discussed[8, 16]. Among those theoretical 
efforts, the possible enhancement from the QCD anomaly is of particular interest since it is 
well known that the rf meson plays an very special role in the dynamics of low energy QCD 



The mass eigenstates rf and rj are related to flavor octet and singlet states rjs, rjo through 
a mixing matrix: 



where 9 is the mixing angle and rjs, ?7o have the flavor content: rjs = -^(uu + dd — 2ss) and 

rj = -^(uu+dd+ss). The associated axial currents are j'5 8 = ■^{u^^^u+d^^^d— 257^755) 

and jg = ^(^7 M 75^ + d'y^jsd + 37^755) respectively. Through the QCD anomaly, the 
divergence of the flavor singlet axial current is non-zero and is given by 



where G^ u = ■^e ,Jvlxr Gp a is the dual of G^ u . This breaks the global U(1)a symmetry for 
massless particles and makes the flavor singlet state t/q evade to be a Goldstone boson of 
chiral SU(3)l <8> SU(3)r symmetry. The QCD anomaly gives its main contribution to the 
large mass of rf (m,< = 0.958GeV) which is much heavier than the other flavor octet states 
such as Ti, Ks and suggests a large gluon content in 7/. 

The QCD anomaly indicates a strong coupling between rf and the gluon field. It is then 
natural to understand the large branching ratio of B — > rfX s through the QCD anomaly. 



[17]. 



rj — r) 8 cos 9 — 7] sin 9, 



r\ = r/ 8 cos 9 + r] sin 9, 
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In the literatures there are two different ways to handle this problem. The one is through 
two-body decay process b — > srj' from some effective Hamiltonian due to QCD anomaly [9]. 
The other one is through three-body process b — > sgrj'[7, 8]. In the first step decay, the 
b quark decays into the s quark and a virtual gluon g*, then g* decays into 7/ and a on 
shell gluon g. This model has some advantages in explaining the spectrum of invariant 
mass distribution of recoiling hadrons. However, the effective gg*rj vertex seems to be too 
small from various approaches [8, 18-20]. In both of the approaches, the effective coupling 
between 7/ and gluon may contain complicate non-perturpative quark-gluon interactions. It 
is then better to treat it as a free phenomenological parameter rather than to evaluate it 
from perturbative calculations [9]. 

In this note we focus on the phenomenological analysis of the first mechanism. The 
effective Lagrangian in this model is given by [9] 

H ef f = aa s G F s L b R G^G^ + h.c (3) 

where a s and Gf are the strong coupling constant and Fermi constant, a is the effective 
coupling. From this effective Hamiltonian, the decay B — > r]'X s arises from the subprocess 
b — > srj'. The evaluation of matrix elements is straightforward: 

(srf\ H eff \b) = aa s G F (s\ s L b R \b) (rf\ G ta/ G tw |0> (4) 

Applying relation (0| j$ 8 ^ \vs{o)) = ifs(o)P >1 to the divergences of both flavor singlet 
and octet axial currents and ignoring small u, d quarks masses, the matrix elements 
(0| G^ U G^ U \rf) can be rewritten as 

(Ol^G^W) = ^y| m ;,(/ 8 sin^ + V2f cos6) 

<0| a s G, v G^ \q) = -^y|mJ(/ 8 cosfl - v^/osinfl) (5) 
In the b quark rest frame, the decay branching ratio is given by 

Br(B^r,'X s ) = ^-T B a 2 G 2 F m*, (f 8 sin 9 + V2f cos O) 2 ^ ~™"' J (6) 



Br(B ^ r]X s ) = ^- TB a 2 G 2 F m' l Jf 8 cos 9 - v^/o sing) 2 ^ ^ (7) 
12 ' mi 

where t b is the lifetime of B meson. 



2. Recoil mass distribution 

For two-body like subprocess such as b — > rj's, the invariant mass is directly related to 
the energy £y of the rj' meson through the relation: m 2 x = m 2 B + m 2 , — 2mg£,', where m B 
and m v i are the masses of B and 7/ meson. The small s quark mass has been ignored. In the 
two-body decay of b — > rj's, the energy of 7/ is fixed from energy-momentum conservation. 
A typical value of the pole mass m& = 4.8 GeV will lead to a narrow peak with the invariant 
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mass of mx s — 1.5 GeV. This seems to be disfavored by the current data since the experiment 
reported a peak at about 2 GeV with a relative wide width in the recoil mass distribution 
[1, 3]. 

However, the above estimation may be too naive. Note that in the two-body decay 
process, the exact distribution of the recoil mass strongly depends on the wave function of 
the B meson which is theoretically hard to estimate. It is too early to draw the conclusion 
that the current data already disfavored all the two-body models. 

To illustrate the non-perturbative bound state effects here we adopt a simple model 
proposed by Altarelli et.al. [21] a number of years ago which is based on the Fermi motion 
of the b quark inside B meson. The basic idea of this model is that the Fermi motion of the 
b quark and the spectator quark q in the B meson make them have back-to-back relative 
three-momenta p in the B rest frame. The momentum is assumed to obey a Gaussian 
distribution as follows 
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0(p) = ^-e-^, p=|p| (8) 



F 



where <p(p) is normalized as J Q °° cf){p)p 2 dp = 1. The mean value of p is < p >= \p F - In this 
model the spectator quark q is on always handled as on shell while the b quark is treated as 
off-shell. Through energy-momentum conservation, the effective mass W of the b quark is 
determined as 



W = m 2 B + m - 2m B \/m 2 + p 2 (9) 



and the energy of the b quark is Ew = \/W 2 + p 2 . Here one parameter p F is introduced 
which specifies both the distribution width and the mean value. As p F is linked to the 
average energy of b quark inside the B meson, in principle it can be calculated from theories 
based on non-perturbative methods or from some models. For example, the calculations 
from QCD sum rule [22] give p F = 0.58 ± 0.06 GeV, and the value from relativistic quark 
model [23] is 0.54 ±0.16 GeV. The value of p F can also be extracted directly from the data. 
A fit to B — > Xg'-f photon energy spectrum give a value of about 0.45 GeV[24] while the fits 
to semi leptonic B decays and B — > J /ipX give a value of 0.57 GeV [25, 26]. Thus the value 
of p F is likely to lie in the range 0.4 < p F < 0.6 GeV. After including the Fermi motion, the 
differential decay width dT(mb)/dmx, should be replaced by 

dT f Pma * , / / \ 2 dY{W) , , 

dp <p{p)p 2 ■ , v ; (10) 



drax Jo dm? 

where p ma x is the allowed maximum value of p, dY{W) / dmx is the differential decay rate 
in the B meson rest frame, which is linked to the one in b quark rest frame through a 
Lorentz boost [25, 27] In Fig.l the invariant mass distribution is generated in this model 
with different values of p F . Here we use ^ which is normalized to unity and independent 
on the value of a. It can be clearly seen that the p F dependence is rather strong. The 
peak of the distributions significantly shifts from around ~ 1.4 GeV ( for p F =0.4 GeV) to 
~ 1.8 GeV (for p F = 0.6 GeV ). Considering the considerable uncertainties in both theory 
and experiment data, there is no significant disagreement in the recoil mass distribution of 
B -> rj'X s . 
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FIG. 1: Recoil mass distribution from process b 
solid, dashed and dotted curves correspond to p F 
at 0.15 GeV. The shadowed area indicates the acceptance cut of mx < 2.35 GeV from the CLEO 
experiment 



0.6, 0.5 and 0.4 GeV. The value of m q is fixed 



3. Bound on a from inclusive and exclusive B decays 

The value of a could be constrained from the exclusive decay modes B — > rj'(rj)K^*\ Note 
that although predictions of the standard effective Hamiltonian approach are too low to 
account for the data of inclusive decay modes, the disagreement in the exclusive decay 
modes are smaller [6, 15]. Furthermore, the effective Hamiltonian approach can reproduce 
the correct patterns oiBr(B —> rj'K) > Br(B —> rjK) and Br(B — > rj'K*) < Br(B —> rjK*) 
which is observed in the experiments. It implies that in exclusive decays modes, it may still 
play an important role, and the interference between different contributions may also be 
important [28, 29]. 

Nevertheless, by saturating the current data of exclusive decays, the upper bound of the 
parameter a can still be obtained. The decay amplitudes of decay modes B — > rfK^ and 
B — > rjK^ in this model read 



M(B 



±,o 



rj'K 



±,0\ 



4tt / 3 

aG F — \l -ml,(f $ sin 6 + a/2/ocosi 



2 2 

m B — m K BK 



2(m b - m s ) 



M{B ±( - 0) ^r/K 



±(0)i 



aG 



4vr 3 



m?(/ 8 cos 9 - v^/osim 



rn 



B 



m 



K 



M(B^ -^r)'K*±W) = - 



2(m b - m s ) 

aG F ^JL% (f 8 sin 9 + V2f cos 9) &^^- A i 
3 V 2 ' (m b + m s ) 



(m b + m s ) 

\P v K*\ m B 



A (mV 



(m b + m s ) v 

where \P v 'k*\ — \PijK*\ — \rfiB- F^ K (q 2 ) and A (q 2 ) are the form factors for B — > K 
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and B — > K* transition with momentum transfer q 2 . The value of mj is taken to be the 
effective one. i.e m^ 2 ~ f W~ 2 (p(p)p 2 dp. In the calculations, we take mj = 4.65 GeV which 
corresponds to p F = 0.5GeV and m q = 0.15GeV. 

The corresponding branching ratio can be evaluated through the following relation 



tb\P\ 



Br =^ M ? (12) 

where \P\ is the momentum of one of the final state mesons in B rest frame. 

It is useful to define two kind of ratios which are independent of the parameter a: 

1) The ratio between B — > rj'X and B — > rjX(X = X s , KorK*). This ratio is independent 

of the value of a and only sensitive to the rj' — rj mixing. In this model one finds [30] 



R 



Br(B r]X s ) Br(B -> rjK) Br(B -> rjK* 



Br(B -> r)'X s ) Br(B -> rj'K) Br(B -> rj'K*) 

m j ( f s cos9- V2f sm9 
mj, I / g sin# + V2f cos6 



(13) 



In the following numerical calculations we take 9 = —17° and fs — fo — 1-06/^ [31] as an 
illustration. This leads to a value of R ~ 0.16. Considering the CLEO data of R < (0.1 ~ 
0.8) [2], it follows that with the constraints from B — > rjK, this model can account for at 
most ~ 60% of the observed B rj'K branching ratio. Note that the exact value of R may 
vary with different sets of parameters 9, f$ and fo] the constraints from R are only an order 
of magnitude estimate. 

2) The ratio between B — > PK* and B — > PK (P = rf or 77). In this model it is 
independent of both the value of a and details of rj' — rj mixing. 

= BrjB -> rj'K*) (m| + m 2 K , - m 2 ) 2 - 4m 2 B m 2 K , 
~ Br(B -> rj'K) (m 2 B - m 2 ^ 2 ' 1 *>bx<~* ^ ' 




The values of F and A in the BSW model [32] are Fq = 0.38, A = 0.32 which corresponds 
to R = 0.84, while from light cone QCD sum rule [33, 34] F = 0.35 ± 0.05, A = 0.39 ± 0.1 
and R — 1.1 ± 0.3. Thus if this model gives the dominant contribution to these modes, the 
value of R should be around 1. However, the current data gives a value of R' < (0.5 ~ 0.4) 
[2]. This is a more clear and stronger constraint than the one from R. With the observed 
small value of R' , this model can explain at most half of the branching ratio of B — > rj'(rj)K 
and therefore is not the dominant mechanism of these processes. In Fig.2(c-f) the numerical 
results of branching ratios as a function of the effective coupling a is given and compared 
with the data. As some inclusive decay modes have not yet been observed by the Babar and 
Belle collaborations, only the CLEO data are used in the numerical evaluations. It can be 
seen from the figure that the data of exclusive decay modes B — > rj'K* and B —>■ rjK impose 
strong constraints on the effective coupling. With these constraints, the maximum value of 
a lies in the range: 

a < (8 ~ 9) x 10" 3 GeV" 1 (15) 

From Eq.(6) and (7), the branching ratio of inclusive decays B —>■ rj'(rj)X s as a function of 
a is plotted in Fig.2(a-b) and compared with the data. In the decay B — > rj'X s the acceptance 
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cut effects is taken into account, which leads to a 19% reduction from the calculation in 
Eq.(6). Given the upper bound of a in Eq.(15) this model can still successfully reproduce 
the B — > f]'X s branching ratio within lcr range. 



4. Prediction of radiative decay J /if) — > 7/777 

From the effective Hamiltonian in Eq.(3), this model can also contribute to the radiative 
J /if) decays into 77'. Using relation Eq.(5), the ratio between J /if) — > 77' 7 and J /if) — > 777 can 
be predicted and found to be the same as in Ref. [31] 



T(J/if) -> 7/7) 
T(J/if) -> 777) 



(0\GG\rj') 



(0\GG\ V ) 



1-m 2 ,) 3 t 

(16) 



which is in good agreement with the data. 

Furthermore, given the value of the effective coupling a the decay rate of J /if) —> rfrpf 
can be predicted. To this end let us first define the ratio 

r{rf) = T(B - rfX 8 )/T(B - g*X s ) (17) 

which can be understood as the size of b — > st]' relative to b — > sg. Taking T(i? — > g*X s ) ~ 
1% and a < 0.008 GeV -1 which comes from the bounds from exclusive decays as an example, 
one finds 

r(n') < 0.045 (18) 

Note that the strong coupling constant in the effective Hamiltonian has been separated 
from the effective coupling a and absorbed in the matrix element of (0| a s GG \rf (77)). It is 
expected that there is no significant running on the value of a from energy scale tub to mji^. 

Since the radiative decay of J /if) — > is dominated by the process J /if) — > g*<7*7, the 
branching ratio of J /if) — > rj'rj'^ can be simply estimated as 

Br(J/if)^ 1 X) VM V ; 

The observation of the process J/if> — > 7X gives Br (J/ if) — > 7X) = (17.0 ± 2.0) x 10 -2 . 
Thus taking r(rf) = 0.045 the maximum branching ratio of J /if) — > v/v/'y is estimated as 

Br(J/if) -> 77V7) ^ 3.4 x 10^ 4 (20) 

The decay rate of J /if) — > v/r] r ) can be estimated as follows 

Br{J/if}^irn) 
Br {J /if) -> 77V7) 

Using the value oi R = 0.16 from Eq.(13) one finds for the maximum branching ratio for 
J/ip —> 7/777 

Br(J/if) -> 77'777) ~ 5.4 x 10" 5 (22) 

Considering the detection efficiency of 77' is about a few percent ( through 7/' — > 7777), 
it may be hard to find a signal of such a decay mode in BES due to limited statistics ( in 
BES 5 x 10 7 J/if> samples are collected). But in the forthcoming CLEO-c project 1 x 10 9 J/^ 
samples are planned to be produced. It is then promising to search for the signal and test 
the predictions from this model in the CLEO-c experiment. 
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FIG. 2: Branching ratios for inclusive and exclusive decay modes as a function of a. 

a) For decay mode B — ► rj X s , the dark and light shadows represent the la (2a) allowed ranges by 
current data. 

b) For decay mode B — > rjX s , the light shadows represent the 90% allowed range. 

c) For decay mode B — > ij'K + , the dark and light shadows represent the 1<t(2<t) allowed ranges. 

d) For decay mode B — > r]K + , the light shadows represent the 90% allowed range. 

e) For decay mode B — > r]'K* + , the light shadows represent the 90% allowed range. 

f) For decay mode B — > r]K* + , the dark and light shadows represent the la (2a) allowed ranges. 



